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ABSTRACT 

X. 

We systematically analyse the gluonic penguin- induced charmless decays B — > Kx<f>{Kx 
denotes the meson state sq (q = u or d) ), in the standard model and the two-Higgs- 
doublet model. These processes, being induced at one-loop level, are of great importance 
in measuring the virtual top quark effect in the standard model, and also in searching for 
the non-standard signals in the low energy region. It is shown that the QCD effect is also 
significant in these processes, as in the weak radiative processes B — > X s 7. We also show 
that the charged Higgs contribution can not provide sizable enhancements for the decays 
B — > Kx4>, in contrast to the decays B — * Kxj- It is also found that processes such as 
B — > .Ki (1400)0 and B — > i^(146O)0 have large branching fractions among B — > Kx4> 
decays. 



^-maiLmasaQauephyas. aichi-edu.ac.jp 



1 Introduction 



Previously rare B decays such as b — > 57, b — > sll, b — > si>V, bq — > II and b — > sg, have 
been the subject of some interest in the literature. The object of such studies has been to 
either allow confirmation of the standard model(SM) or to find indications of additional 
effects beyond the SM. Taking such an approach in our previous paper [jj (referred to as 
paper I), we analysed the inclusive decays B — > X s ^ and exclusive decays B —* Kxj 
(where Kx denotes the meson state sq(q=uord) ). In this work we included the 
nonstandard physical effects due to the charged Higgs contribution in the two-Higgs- 
doublet model(THDM). We showed that sizable enhancements for these processes are 
possible compared to the SM. It is well known that, in these processes QCD corrections 
play an important role |2j in the quantitative discussion, resulting in an enhancement of 
three to four times in the branching ratio for b — > 57. 

In the present paper, we discuss another interesting class of processes: B — > K x <t>- 
There is no previous work with full QCD corrections and including higher .ff -meson res- 
onances. The characteristic feature of these processes is that they are mainly induced 
through the gluonic penguin interaction via b — > s g*, where g* denotes the virtually 
emitted gluon. The non-leptonic 5-decay processes to the final states including s + s + s 
have previously been discussed by three of the present authors and others, with the aim 
of clarifying the nonstandard physical effects due to the charged Higgs contribution in the 
THDM [g]. However, QCD corrections were not fully included in these works, and the 
decays into higher X-resonances with the 0-meson were not included. 

Now we want to make a more detailed analysis of the processes B — > Kx<t>- First we 
summarize the QCD corrections to the effective Hamiltonian for the process b — > s + s + s. 
For this process, the QCD corrections induced by the dominant contribution of the t-quark 
up to the one-loop level essentially require the additional dimension-six local operators of 
Wilson's operator-product expansion. These operators have been approximately neglected 
in the analyses of the other rare B decays ||. Such consideration of more complete QCD 
corrections has already been given by Buchalla, Buras and Harlander 0, who provided 
a detailed renormalization group analysis to investigate the ratio e'/e systematically by 
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including the contributions from the gluonic, photonic, Z°, and neutral Higgs penguin 
diagrams. They obtained the Wilson coefficient functions with full 0(otQED) contributions 
in a compact and transparent form. 

The concrete and detailed description of our evaluation of the relevant Wilson functions 
of the linearly independent dimension six operators is given in section 2, although our 
evaluation of them is partly the same as that made by Buchalla, Buras and Harlander. 
Note though that the estimation of the Wilson coefficient functions for full operators in 
THDM has not been previously given. In section 3, we describe the calculation of the 
amplitudes and the decay rates for B — > Kx4>- We emphasize that our calculation of the 
amplitudes uses the specific form factors given by Isgur, Scora, Grinstein and Wise ||, as 
well as the factorization assumption as seen in this section. Our results are compared with 
another model of form factors given by Bauer, Stech and Wirbel |7j . The effect of possible 
charged Higgs contributions is also studied. Recently, new experimental results, in the 
form of a limit on and a non-zero value of respectively, for B(b —>■ sj) and B(B — > K*j), 
have been reported by the CLEO experiment ||. These results lead to new constraints 
on the values of the relevant parameters occurring in expressions for the amplitudes of 
these decays. Section 4 is devoted to brief comments on the resultant constraints. In 
particular, we discuss the constraint on cot/5 versus win in the THDM. In section 5, our 
summary and conclusions are presented. 

2 QCD effects in the rare decay processes of B-meson 

The rare decay process of the 6-quark to the final state s + s + s is induced at the one-loop 
level mainly through b s + g*. 

The fundamental four quark interaction mediated by the process b — > s + g*, i.e. 
induced by the gluonic penguin diagram, is described by the effective Hamiltonian 

where 

rf = G 1 {x t ) lll {\ - 75 ) + 3^G 2 (x t K^(l + 75) (2) 
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and 



rf = F l {y) llx {l - 75) + 3i^F 2 (y)a^(l + 75 ). (3) 



The functions Gi and Fi are given by 



Gi(x t ) 



x t (l - x t )(18 - llx t - x 2 ) - 2(4 - \Qx t + 9a£) lnx 4 



;i - ^) 4 

,2 „2j 



r f l , , \ rnl- qH(l -t) 5 , 

+ 48{ / dt t(l - t) In — \ - : '- r } 

Vo V ; m 2 w (l - t) + m 2 c t - q 2 {\ - t) 36 J 



(1 - x t )(2 + hx t - x 2 ) + 6x t lnx t , A , 
G 2 (x t ) = x t , (4) 

F M = ycot ^(_^yM_-29y + 7y^ +6 (2-3y)lny 



3(1-2/) 
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^2(2/) = -cot (3 G 2 (y) + 2y — — , 

where x 4 = m^/m^, y = m 2 H /m 2 and cot /3 = Urf/v u is the usual notation for the ratio of 
the vacuum expectation values of the neutral sectors of the two-Higgs-doublets. Here, in 
the first equation for Gi(x t ), we assume V t b = V cs , V ts = V c b and V u b = . Note that we 
have chosen a specific variant of the 2HDM, i.e. of the form of the coupling of the Higgs 
doublets to the fermions. There are four variants of the 2HDM , distinguished by different 
schemes for the Yukawa couplings to the fermion sector. The reader will find a concise 



summary of the possible models in reference |TIJ. For what follows we will take model II 



of that work, in which one doublet gives mass to the up-type quarks, and the other to the 
down-type quarks and charged leptons. As we will discuss, there are considerable uncer- 
tainties in the calculations that follow, and as a consequence, we will content ourselves to 
investigate the possibility of distinguishing between SM results and those with contribu- 
tions arising from charged scalars. Distinguishing further between competing models will 
in general require considerable refinements in the theoretical calculations. We also note 
that, in the second term of the r.h.s., we should include the charm quark contribution due 
to the soft GIM cancellation. This term reduces to the familiar result 8 In m 2 ./m 2 v in the 
limit of the squared momentum of virtual gluon vanishing; q 2 = 0. Note though that this 
limit is not realistic for the virtual gluonic penguin induced processes, since the gluon is 
off mass-shell. In the analysis to follow we take the appropriate nonzero value of q 2 . 



The following four-quark operators and the magnetic-transition-type operators are 
relevant for the processes under consideration. We write the Hamiltonian 

Heff = ^VttKECiMOitjJ.) , (5) 

where \i in the parentheses denotes the energy scale at which the operators are relevant 
for the decays. The operators Ot are defined as follows: 

3 = (s La 'fb La )( QlpI/hQlp) , 

5 quarks 

C 4 = (s L aYb L p){ VLfSl^La) , (6) 

5 quarks 

5 quarks 

5 quarks 

g 

°7 = -i-^ m bSL .o llv b Ra (i ll e v , 
°s = -t^m b s La a^T^b R ^e a u . 

Coefficients relevant to our processes are defined at the energy scale m w to be 

Ci(m w ) = , 

C 2 (m w ) = 1 , 

1 1 
C 3 (m w ) = --C A (m w ) = C 5 (m w ) = --C 6 (m w ) 

° ,( "" v) (C lW + F lW ) 1 (7) 



288vr 

C 7 (m w ) = C 7 SM (x t )+C 2 7 H (y) , 
C 8 (m w ) = - l (G 2 {x t )+F 2 {y)). 



Here Cf M (x t ) and C^ H (y) are 
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and 

C 2 7 H (y) = 77 V V J %2 ~ % + 3 --(3^-2)1^] 

a„ y r 8y 3 - 3y 2 - 12y + 7 2 . 3 . 

respectively. In Eqs.(7), we neglect the contributions from Z°- and photonic-penguin 
interactions [O] since these terms are very small numerically due to a s < a and do not 
affect the present analyses. Also the second term, i.e. the charm contribution, in G\{xt) 
should be taken into account |L2| at the scale \i = m b by setting q 2 ~ ml/2 JK 



Numerically, in the energy region between m\y and mj, we put the flavour number 
/ = 5 in the renormalization group equation |f[4 |. We evolve the coefficients Cj(//) starting 



from the scale mw as given in Eqs.(7) to the scale \i = m b = 4.58GeV, and then we obtain 
the following numerical coefficients: 

Ci(m 6 ) = -0.240 , C 2 (m b ) = 1.103 , 

C 3 (m b ) = 0.011 + 1.125C 3 {m w ) -0.121C 4 (m w ) , 

C 4 (m 6 ) = -0.025 - 0.291C 3 (m w ) + 0.82AC A (m w ) , 

C 5 (m b ) = 0.007 + 0.9UC 3 {m w ) + 0.083C 4 (m w ) , (10) 

C 6 (m b ) = -0.030 + 0.229C 3 (m w ) + 1.465C 4 (m w ) , 

C 7 (m b ) = -0.199 + 0.629C 3 (m w ) + 0.931C 4 (m w ) + 0.675C 7 (m w ) 

+ 0.091C 8 (m w ) , 

C 8 (m b ) = -0.096 - 0.598C 3 (m w ) + 1.029C 4 (raw) + 0.709C 8 (mvK). 

In Table 1 we summarize the numerical values of the above coefficients in the THDM for 
two values of the parameters (ma, cot/3) [denoted in table by 2HDM(1) and 2HDM(2)], 
as well as the SM results for typical values of the SM parameters V^,- and m t . In order to 
calculate the enhancement due to QCD corrections, we have compared the coefficients in 
both cases with and without these corrections. Note that most of the previous predictions 
for B — > Kef) and K*(p decays have been given without the QCD correction. 
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Table 1 



In Table 1, the heading "Without" QCD means that we include no RG evolution of the 
coefficients and fix their values at the scale raw, as given in Eq.(7), with q 2 ~ ml/2. We 
see that the values of Ci(rrib) with the QCD correction are larger than those without by 
a factor of 1.5 ~ 2. The values Ci(mfe) ~ C^m^) are almost identical in the SM, the 
THDM(l) and the THDM(2). It is noteworthy that the differences between the three 
models exist only for the coefficients CV(mb) and ^(mi,). As mentioned later in section 
4, the value of cot (3 is strongly limited by the recent CLEO experimental search || for 
the B — > Xg'-f decay. As a result, here we use the values of m# and cot (3 which satisfy 
the present experimental restrictions. In the next section we apply the above analysis on 
QCD corrections to the calculation of the decay amplitudes. Note that no clear differences 
among the SM and the THDMs for the coefficients C\{m^) ~ C^m^) are seen in Table 
1 even using the QCD corrected coefficients. So we expect that the decay B — > X s 7 will 
be the best process to search for an effect beyond the SM. However, it is still important 
to analyse the gluonic penguin effect through the exclusive processes induced by this 
interaction. This allows us to check the one-loop effects based on the SM, or to study the 
signals from beyond SM physics. Such an analysis is given in section 3. 

3 Exclusive decays B — > Kx4> induced by the gluonic 
penguin interaction 

Of the penguin induced charmless B meson decays, the exclusive non-leptonic rare 
decays B — > K x <t> are of special interest, since these processes are the typical of those 
caused by the gluonic penguin interaction. Analysis of these exclusive decays is also helpful 
for the future experimental study at 5-factories. In this section, we calculate these decay 
rates by including QCD corrections described in section 2. The contribution from new 
physics is also studied. Since the form factors play an important role in estimating the 
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branching ratios precisely, we also study the form factor dependence of our results. 

In the sq system, a rich spectrum of states has been observed. The resonance state 
Kx is specified by the quantum numbers n, L, s and J, which denote the radial excitation 
quantum number, the orbital angular momentum, the sum of the spins of the two quarks 
and the total spin of the meson, respectively. We investigate the following eight states 
with the notation n 2S+1 Lj: 1%, 1 3 S U 1 3 P 2 , l 3 Pi, l x Pi , 1 3 P , 2 1 S and 2 3 S 1 . Physical 
mesons corresponding to each state are p 



l 1 ^ : K , l 3 ^ : iT(892) , 1 3 P 2 : K*(U30) , 

l 3 Pi ,1^1 : Ki(l270) , Ki(1400) , (11) 

1 3 P : K (U30) , 2 1 5'o : #(1460) , 2 3 S 1 : AT* (1410) . 

The spin 1 mesons i^i(1270) and Ai(1400) are nearly 45° mixed states of l 3 Pi and l 1 Pi- 

We proceed by showing the formulation to calculate the decay matrix elements. The 
operators 3 , 4 , 5 , 6 and 0' 8 are relevant for the decays B — > K x 4>- Here the operator 
Og is derived from the interaction in which the virtual gluon in the magnetic operator Og 
couples with s + s, to give 



a 1 

O's = -im b ^-s La a^T^b Rp q,Sp lv T« a s a . (12) 



By the use of these operators, the decay amplitude may be written as 

AC 1 

{K X (f) | H e}f I B) = -y=-V tb V* £ Ci(n)(Kx4> | Oiifi) | B) . (13) 

3,4,5,6,8' 

In the following analysis we take the value of the coefficient C 8 (yu) as equal to Cs(/x). 
Although Cs(/i) does not include full the QCD correction of Cg(/i) in the leading log 
approximation, this replacement does not seriously affect the results numerically, as the 
0' 8 term is the next to leading one, compared to the operators O3, O4, O5 and CV 

We also use the factorization approximation in order to estimate the hadronic matrix 
element. The factorization assumption works successfully in D meson and B meson decays 
TT| within a factor two in the branching ratios. Then, the hadronic matrix elements of 
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the above operators are given by 

(K x <f)\0 3 \B} = (K x( f)\0 4 \B) = ~(K x <j ) \0 5 \B)=4(K x <j)\0 6 \B) 

= ^((J>\s^s\0)(K x \sY(l-l 5 )b\B}. (14) 

Using the Gordon identity and the color algebra relation, the operator 0' 8 reduces to 

a 1 8 4 
0' 8 = i^m b - 2 [i-m b {s Ll »b L ){s Ll »s L ) - -(sa, u b)s(p') L W - tW] (15) 

in the limit m s = 0, where a Fierz transformation is performed and the relation (0 | slSr \ 
0) = is used. Then, the hadronic matrix element is factorized as follows: 



(K x( p \0' S \B) = i^ mb ±-{£- mb (<p | s Ll Sl | 0)(K X \ s L ^b L \ B) 

47T q z 9 
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with 



\(<t> | s l (Yp' u ~ 7V)si I Q)(K X I sa^b \ B)} , (16) 



(0I*l(tV*-7VV|o) = (17) 

and where g^ = 0.23GeV 2 is taken and p /fl = p^/2 is assumed. 

The hadronic matrix elements (K x \ sa^b \ B) are given for the l 1 ^ and l 3 ^! states 
in terms of the form factors by 

(K(p') | sa^b | B{p)) = is T [{p + p%(p - p') v - (p - p%{p + p') u ] , 
(K*(p', e) | sa^b | B(p)) = gle, uXa e*\p + p'Y + g%vx*e*\p - p'Y (18) 

+ h T e, uXa (p + p')\p- P r^-p) . 

The form factors appearing in the above equations are easily estimated by using the 



relations derived from heavy quark symmetry [17 

T U ' f- 



S 



2rrih 



m b 2m b 

9 T + -g T = -2m b g, (19) 

T I T f . P-P' 

9+ + 9- = — + 2 g, 

m b m b 



where the form factors /+, /_, g, a + and a_ are defined in the hadronic matrix elements 
of the vector and the axial- vector current JT7j. For the 1 3 P2, 2 1 S , and 2 3 S' 1 states, similar 
relations are satisfied, while for the l 3 Pi, l l P\ and 1 3 Pq states, the overall signs of r.h.s. 
in Eq.(19) should be reversed. 

These form factors are to be evaluated at (p — p') 2 = m 2 ^ for each final state and they 
generally depend on the quark potential. Here we use the form factors given by Isgur, 
Scora, Grinstein and Wise || , which have been successfully applied to the electron energy 
spectra of the semileptonic D and B meson decays. These are derived using harmonic 
oscillator type wave functions with a variational method, assuming a Coulomb plus linear 
potential. This simple model gives quite reasonable spin-averaged spectra of cd and bd 
mesons up to L = 2. These form factors include the relativistic compensation factor, 
although their model itself is a nonrelativistic one. 

The decay amplitudes are given by 



AG 

(K x <p \ H eff \ B) = -j=-V tb V*{A^(K x | s Llll b L | B) 

+ iB( V %-p^){K x | sa^b | B)} , (20) 



with 



2 ix q z 



B EE ^Cs(fl), (21) 
47r q z 

and where p is the three momentum of the K x meson. 

The decay widths of the possible modes are given in terms of the form factors by 



V(B K x <>) | V tb V* | 2 p( A' 2 X A - B 2 Xjj ! ADX. W ). (22) 

A7xm B 



The quantities X^, Xb and Xab for each state are given in terms of the form factors as 
follows; 

for the 0~(l 1 S'o) state, 

X A = ^/ip 2 , X B = WmlmyV, X AB = -8m 2 B f + s T p 2 , (23) 
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for the l"(l 3 5'i) state, 

X A = \{f + 4m!sV ) + } 2 {(p K * ■ P4> )f + 2m 2 B a + p 2 } 2 , 

X B = 32m|^V, X A B = 8m 2 B gglp 2 , (24) 

for the 2+(l 3 P 2 ) state, 

XA = ^-(k 2 + AmlhV)p 2 + -^{(p K2 .p,)k + 2rn%b +P 2 }Y, 
^m K% om K 2 m <f> 

4 4 
m B J2J v /( m B i T „4 



X B = X AB = 4-fHV. (25) 

for the l + (l 3 Pi) or the l+^Pi) states, 



m K 2 m K 2 



X A = \{t 2 + 4m|gV) + 1 2 {(p Kl ■ Plf) )i + 2m 2 B c + p 2 } 2 , 
2 4m z Kl m% 



X B = 8{m^gi + -gi_)-2(p Kl -p^gi + } 2 

+ — 2 i{ m i(9i+ ~ 9i-) ~ 2 (?Ki ■ P*)0+i}(p*i ■ P<t>) 



+ 2m 2 B {g T 1+ - mlhl)p 2 } 2 , 



2 

XaB = o 2"{(P^i ■ P4>)^ + 2m B c +P 2 } X 

[{(fl , i+ ~ 9i~)m 2 - 2(p Kl ■ p^)gi + }(p Kl ■ P^>) + 2m 2 B (gf + - m 2 h^)p 2 ] 
- 4£{ml(gT + -gl)-2(p Kl -p (f> )gl + }, (26) 

and for the + (l 3 P ) state, 



m 



X A = ^u 2 + p 2 , X B = 0, X AB = 0. (27) 

For the _ (2 1 S , ) and l~(2 3 5'i) states, similar forms of Xa, X b and Xa B as for _ (1 1 S , ) 
and l~(l 3 Si) apply, respectively. The explicit forms of the form factors /+, /, g, a + , k, 
h, b + , £, q, c + and u + are given in the Appendix of Ref.0. The form factors s T , g+, gj + , 
g% + and hj are calculated by using Eqs.(19) and the form factors referred to above. 

The decay branching ratios are shown in Table 2 in the SM, where Cj(/i) are evaluated 
at /J, — mi, — 4.58GeV. 



Table 2 
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It is found that the decays B — ► K\ (1400)0 and B — > ff(146O)0 dominate the total decay 
rate of B — > Kx4>- We also show the predicted values "without" QCD corrections, which 
cases were investigated by two of the present authors previously in the second reference of 
Refs. 0. It should be remarked that the QCD effects significantly increase the branching 
ratios by almost a factor of three. This is to be compared to the process b — > sj, in which 
the QCD corrections enhance the decay ratio by one order as discussed in paper I. Since 
the predicted values depend on the values of the form factors, we compare our predicted 
branching ratios for K<p and K* (890)0 final states with the ones calculated using another 
model of the form factors given by Bauer, Stech and Wirbel J7|. As shown in Table 3, the 
predicted values in the latter model are almost three times larger than the ones in the 
GSIW model. Thus, the predicted values depend significantly on the model of the form 
factors, and it would obviously be of great value to improve the reliability of choosing 
these form factors. 



Table 3 



These form factors will possibly be tested by the radiative rare decay B — > K* (890)7, 
which is the process discussed in the next section. 

Finally, we study the effect of new physics on these decays. As a typical example of 
such new physics, the effect of charged Higgs bosons is investigated. Then, Ci(/z) are 
modified as discussed in section 2. We show the predicted branching ratios for the typical 
case where m# = lOOGeV and cot /3 = 1 in Table 2, where we see that the effect of the 
charged Higgs boson increases the decay rate at most by 10 ~ 18%. This is the extreme 
case to see the Higgs effect and seems to be excluded by recent CLEO experiment [|J as 
discussed in the next section. Then the charged Higgs boson appears to have a minor 
contribution to these decays. Thus, it may be difficult to observe the evidence of the new 
physics in these decays if we take into account the large ambiguity in the choice of the form 
factors, as well as the assumption of factorization. However, the radiative exclusive rare 
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decays such as B — > K* (890)7 may help to choose between models of the form factors, and 
the factorization assumption will be studied precisely in the normal non-leptonic decays 
of the B meson. Then the theoretical calculations of these rare decays will become more 
reliable, and the contribution of effects due to new physics will be able to be discussed 
quantitatively. 



4 The renewed analysis based on the new CLEO ex- 
perimental results 

Recently, the new experimental limit on B(b — *> 57) and the new value of B{B — > K*^) 
have been reported by the CLEO group as follows; 

B(B -> X s7 ) < 5.4 x 10~ 4 , (28) 

B{B -> K*i) = (4.5 ± 1.5 ± 0.9) x 10" 5 . (29) 



The experimental data on these processes improve the previous constraints |18| on the 
charged Higgs mass and the parameter cot (5 in the THDM. In Fig.l, we show the modified 
allowed parameter region of cot (3 versus m# obtained from the new upper limit of B(B — > 
Xgj), which replaces that of the Fig. 2 of paper I. 



Fig.l 



The value of m# is restricted to be larger than about 220Gev which is for cot (3 = 0, m t = 
140GeV and = 4.58GeV. It is of interest to note that the new experimental upper 
bound B(B —>■ X s ^) < 5.4 x 10~ 4 is very close to the QCD corrected predicted value of 
3.3 x 10~ 4 (for m t = 140GeV) in the SM. 
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The predicted branching fractions of B — > K* (892)7 given in the paper I are 
' 3.0 x 10" 5 



B(B -> AT* (892)7) 



(5.3 ~ 6.9) x KT 5 (TBDM:m H = 200GeV, cot (3 = 0.0 ~ 2.0) 

(30) 



at m t = 140GeV. These values are consistent with the recent CLEO experimental results 
given in Eq.(29). At present it is hard to distinguish the effect of the charged Higgs boson, 
because of large experimental errors. In our calculation, for both cases, the ratio of the 
exclusive decay rate of B — > K* (892)7 to the inclusive decay rate of B — > X s ^ is 7%. 

Further experimental progress on both the inclusive and the exclusive processes is 
eagerly anticipated, to help settle whether the SM correctly describes these decays, or 
whether we should take into account the possibility of physics beyond the SM. 



5 Summary and conclusion 

We have analysed the decays B — > Kx4> i n the SM, and also studied the effect of 
inclusion of the charged Higgs contribution. 

First we presented the QCD effects for the effective Hamiltonian relevant for b — > 
s + s + s, which is mainly induced at the one-loop level through the gluonic penguin, 
where the approximately neglected operators in the previous analyses of the other rare B 
decays are properly included. Then the numerical values of the relevant Wilson coefficient 
functions Ci(rrib) in the SM as well as in the THDM are summarized, where we find that 
the noteworthy differences among the SM and the THDMs are seen only for the two 
coefficients CV(mb) and Cs(m&). 

In estimating the hadronic matrix elements of the exclusive decays B — > Kx4>, fac- 
torization and a specific model for the form factors are assumed, the latter being those 
given in Ref. 0, based on the harmonic oscillator type wave functions with variational 
method using the Coulomb plus linear potential. As well, these form factors are based on 
relations derived from heavy quark symmetry. 

We obtained the branching ratios of these processes in the standard model and then 
we found that the decays B — > K\ (1400)0 and B — > if (1460)0 are the dominant modes 
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of B — > Kx4>- We also found that, as in the weak radiative processes B — > -Xyy, the QCD 
effect significantly increases the branching ratios by a factor of about three. However, 
the predicted values depend significantly on the model of the form factors, which we have 
shown by comparing the results obtained by using two sets of form factors, those of Ref . || 
and Ref. 0. Accurate measurements of decays such as B — > K* (890)7 ma y be helpful to 
select the most reliable model of the form factors. We also studied the effect of the possible 
existence of charged Higgs boson contributions to these decays. However, according to 
our calculation, it gives only a minor contribution, increasing the decay rates at most 
by 10-18%. Due to the theoretical ambiguities of the prediction, such as its dependence 
on the model of form factors and the factorization assumption, such minor contribution 
makes the clear observation of evidence for the existence of the charged Higgs contribution 
difficult. So weak radiative decays may be more favourable to search for charged Higgs 
effects than the processes we consider here, B — > K x 4>- 

note added: After most of this work was completed, we have received the related 



paper by R. Fleischer [12 
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Table Captions 

Table 1: Co efficients C^rrib) to Cg(mb) of the interactions Eqs.(6) in the SM and 
THDM "without" and "with" the QCD effect. C x (m h ) and C 2 {m b ) are given in Eq.(lO). 
The mass of top quark is fixed to be 140Gev and mfc=4.58GeV. The THDM(l) and 
THDM (2) results correspond to the case m# =220GeV, cot (3 — and txih = 300GeV, 
cot/3 = 1, respectively. 

Table 2:Branching ratios in the standard model "with" and "without" QCD ef- 
fect. The predictions including the contribution from the charged Higgs bosons (m# = 
100GeV,cot (3 = 1) "with" QCD are also shown in the last column. 

Table 3:Predicted branching ratios using the ISGW form factor model || and the 
BSW form factor model [0 in the SM including QCD corrections. 

Figure Captions 

Fig. 1 The allowed region of parameters m# and cot (3 in the THDM. The allowed 
region is to the right hand side of the corresponding line. The thick and the thin lines 
correspond to the cases m& = 4.58GeV and m b = 5.00GeV, respectively. The three cases 
m t = HOGeV, 140GeV and 170GeV are shown as indicated in the figure. 
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